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This paper is the fourth in a series devoted to identifying and explaining the properties of strongly 
correlating liquids, i.e., liquids where virial and potential energy correlate better than 90% in their 
thermal equilibrium fluctuations in the NVT ensemble. For such liquids we here introduce the 
concept of "isomorphic" curves in the state diagram. A number of thermodynamic, static, and 
dynamic isomorph invariants are identified. These include the excess entropy, the isochoric specific 
heat, reduced-unit static and dynamic correlation functions, as well as reduced-unit transport coeffi- 
cients. The dynamic invariants apply for both Newtonian and Brownian dynamics. It is shown that 
after a jump between isomorphic state points the system is instantaneously in thermal equilibrium; 
consequences of this for generic aging experiments are discussed. Selected isomorph predictions are 
validated by computer simulations of the Kob-Andersen binary Lennard- Jones mixture, which is 
a strongly correlating liquid. The final section of the paper relates the isomorph concept to phe- 
nomenological melting rules, Rosenfeld's excess entropy scaling, Young and Andersen's approximate 
scaling principle, and the two-order parameter maps of Debenedetti and coworkers. This section 
also shows how the existence of isomorphs implies an "isomorph filter" for theories for the non- 
Arrhenius temperature dependence of viscous liquids' relaxation time, as well as explains isochronal 
superposition for strongly correlating viscous liquids. 



How much does knowledge of a system's thermal equilibrium fluctuations at one state point tell us about its 
behavior at other state points? In principle, complete knowledge of the fluctuations provides enough information to 
determine the density of states, from which the free energy at other state points may be calculated. In practice, 
only second-order moments of the fluctuations may be determined reliably These generally give little knowledge of 
the system's properties away from the state point in question. It was recently shown that a large class of liquids 
exhibit strong correlations between their virial and potential energy (NVT) thermal equilibrium fluctuations.* 1 * 2 * 3 * 4 ^ 5 -*^ 
Such liquids have a hidden (approximate) scale invariance.* 5 -*-^ Because of this, important global information about 
the system may be obtained from knowledge of the virial and potential energy fluctuations' second-order moments at 
one state point. This unusual situation in statistical mechanics is the background of the present paper, which is the 
fourth in a scries^* 5 -* devoted to identifying and explaining the properties of strongly correlating liquids. 

Paper P^of the series presented results from computer simulations of 13 different liquids. The results show that van 
der Waals and metallic liquids are strongly correlating, whereas hydrogen-bonding liquids like methanol and water are 
not. Likewise, covalent and ionic liquids are not expected to be strongly correlating because competing interactions 
generally spoil the correlations. Paper II 4 gave a thorough analysis of the cause of the strong correlations, which 
is briefly recapitulated below. It was shown how to qualify the simple explanation of the correlations given in our 
first publication (Ref. [TJ, where strong correlations were argued to derive from particle-particle close encounters, thus 
probing only the repulsive part of the potential which is in many cases well approximated by an inverse power law. This 
explanation must be qualified in order to explain the occurrence of strong correlations at low temperatures and/or low 
pressures, as well as in the crystalline state. A number of consequences of strong virial / potential energy correlations 
were also discussed in Paper II. Paper III, 5 published in tandem with this paper, gives further theoretical results on 
the statistical mechanics and thermodynamics of the hidden scale invariance that characterizes strongly correlating 
liquids. Paper III also presents results from computer simulations demonstrating that strong virial-potential energy 
correlations are present even in non-equilibrium situations. 

The present paper introduces the new concept of "isomorphs" in the state diagram of a strongly correlating liquid 
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and derives a number of isomorph characteristics. The existence of isomorphs distils the properties of strongly 
correlating liquids into one single concept and its immediate consequences, a concept that is defined without reference 
to correlations. 

In order to recapitulate the definition of a strongly correlating liquid, recall 8 -* that for a system of N particles 
in volume V at temperature T, the pressure p is a sum of the ideal gas term NksT/V and a term reflecting the 
interactions, W/V, where W is the so-called virial, i.e., 



P V = Nk B T + W. (1) 

This equation is usually thought of as describing thermodynamic averages, but it also applies for the instantaneous 
values. The instantaneous ideal-gas pressure term is a function of the particle momenta, giving the NksT term. The 
instantaneous virial W is the function of particle positions definecPby W(ri, ■■■,r N ) = — l/3J2i r i ' V rj i7(rx, ..., r jv) 
where U(r±, rjv) is the potential energy function. If AU is the instantaneous potential energy minus its thermody- 
namic average and AW the same for the virial, at any given state point the WU correlation coefficient R is defined 
by (sharp brackets here and henceforth denote equilibrium NVT ensemble averages) 



R = ( AWAU ) (2) 
^/((AW) 2 )((AU) 2 ) • 

By the Cauchy-Schwarz inequality the correlation coefficient obeys — 1 < R < 1. We define strongly correlating liquids 
by the condition R > 0.9P The correlation coefficient is state-point dependent, but for all liquids we have studied 
by computer simulatiorP R is either above 0.9 in a large part of the state diagram, or not at all. In all cases the 
correlation coefficient quickly decreases when pressure becomes negative, but strongly correlating liquids generally 
remain so at zero pressure. 

Strongly correlating liquids include^ 2 * 3 * 4 * 6 * 9 * the standard Lennard- Jones (LJ) liquid, the Kob-Andersen binary 
LJ mixture as well as other binary LJ-type mixtures, a "dumbbell" liquid of two different LJ spheres with fixed 
bond length, a system with exponential repulsion, a seven-site united-atom toluene model, the Lewis- Wahnstrom 
orthoterphenyl model, and an attractive square-well binary mixture. Strongly correlating liqui ds have simpler physics 
than liquids in general. This fact has particular significance for the highly viscous phase H | i2 | i3 | i4 | i5 | i6 | i7 | i8 | i9 | 20 | 2i ] 
Thus to a good approximation stro ngly cor relating viscous liquids have all eight fundamental frequency-dependent 
thermoviscoelastic response function d 22 * 23 * 24 * given in terms of just oneP i.e, they are almost single-parameter liquids in 
the sense of having dynamic Prigogine-Defay raticE^l close to un itjl214123| Moreover, strongly correlating viscous liquids 
obey density scaling, i.e., t heir average relaxation time r varies with density p — N/V and temperature according to 
t = F(p 1 /T) j 7 l 25 l 26 l 27 l 28 l 29 l Even complex systems like biomembranes may exhibit significant correlations for their slow 
thermodynamic degrees of freedom; this was shown by all-atom computer simulations of five phospholipid membranes 
which exhibit strong correlations of the energy-volume fluctuations in the NpT ensembleP^ 

When instantaneous values of virial and potential energy are plotted against e ach oth er for a strongly correlating 
li quid in thermal eq uilibrium at constant volume, an elongated ellipse appears? 1 * 3 * 4 * 6 ! The slope of this ellipse is 
y/ ((AW) 2 ) I ((AU) 2 ) . As detailed in Sec. IID this quantity, which is weakly state-point dependent, is to a good 
approximation the exponent 7 of the density-scaling relation, t — F(p~ 1 /J 1 ).* 7 -* 9 -* Thus for a strongly correlating liquid 
knowledge of the equilibrium fluctuations at one state point provides information about how the relaxation time varies 
with density and temperature. 

What causes strong t he WU correlations of some liquids? A hint comes from the well-known fact that an inverse 
power-law pair potential, 1 31 * 32 * 33 * 34 * 35 * 36 * 37 * 38 * 39 * 40 * 41 *^ f(r) cx r~ n where r is the distance between two particles, implies 
100% correlation.™ In this case the slope is rt/3.* 11 ^ In simulations of the standard LJ liquid we found slopes around 
6, corresponding ton= 18.^ Although this may seem puzzling given the expression defining the LJ potential vlj(t) = 
4e[(r /<t)~ 12 — (r/er)~ 6 ], if on e wish es to fit the repulsive part of the LJ potential by an inverse power law, an exponent 
around 18 is indeed required! 1 * 4 * 43 * The reason is that the attractive r~ e term makes the repulsion considerably steeper 
than the bare repulsive r -12 term would imply. 

Paper II gave a thorough discussion of the WU correlations with a focus on the standard single-component LJ 
liquid; this included also a treatment of the (classical) LJ crystal where one finds 0.99 < R < 1 at low temperature. 
According to Paper II the r-dependent effective exponent n which controls the correlation is not simply that coming 
from fitting the repulsive part of the potential, but rather n^ 2 \r) = —2 — rv"'(r)/v"(r). This number is 18-19 around 
the LJ minimum. In fact, the LJ potential may here be fitted very well with an "extended inverse power-law" (eIPL) 
potential)^ a potential of the form v^(r) = Ar~ n + B + Cr with n of order 18. For this potential n^ 2 \r) = n. At 
constant volume the linear term contributes little to the virial and potential-energy fluctuations: When one nearest- 
neighbor interatomic distance increases, another decreases in such a way that their sum remains almost constantPThis 
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means that virtually correct canonical probabilities are arrived at by using the inverse power-law (IPL) approximation, 
an observation which inspired us to the below isomorph definition. 

For an IPL liquid several quantities are invariant along the curves in the phase diagram given by p n /' 3 /T = Const. 
Paper III^I summarizes the thermodynamic IPL invariants, which include Hclmholtz free energy over temperature, 
excess entropy, average potential energy over temperature, isothermal bulk modulus over density times temperature, 
and virial over temperature. In reduced units the dynamics of an IPL liquid is also invariant along the p"/ 3 /T = Const, 
curves!^ The present paper shows that some IPL invariants give rise to general "isomorph invariants" of strongly 
correlating liquids. Not all IPL invariants generalize, however, and e.g. the equation of state of a strongly correlating 
liquid is usually poorly represented by the IPL approximation.^ 

We demonstrate below several implications of one single assumption: the existence of curves in the phase diagram on 
which for any two state points there is a one-to-one correspondence between their respective microscopic configurations, 
such that corresponding configurations have identical configurational NVT canonical probabilities. These curves in the 
state diagrams are referred to as isomorphs. Section II defines isomorphs and summarizes their properties classified 
into thermodynamic, structural, equilibrium dynamic, and aging properties. Most isomorph properties come in the 
form of isomorph invariants. In Sec. II we also discuss how to identify isomorphic curves in the state diagram. Section 
III presents results from computer simulations of (mainly) the Kob- Andersen binary LJ mixture, validating some of 
the isomorph predictions. Section IV relates the isomorph concept to selected topics of current liquid state theory 
and experiment. Section V gives a brief conclusion. 

II. ISOMORPHS 

This section introduces the concept of isomorphs in the state diagram of a strongly correlating liquid. Although the 
definition of an isomorph refers neither to IPL potentials nor to strong WU correlations, only strongly correlating liq- 
uids have isomorphs. This is bec ause the existence of isomorphs reflects the hidden scale invariance that characterizes 
strongly correlating liquids 

A. Isomorph definition 

Assuming that the origin of the coordinate system is centred in the liquid, for any microscopic configuration 
(i"i, ... , rjv) of a thermodynamic state point with density p, the "reduced" (dimensionless) coordinates are defined by 

r, = p 1/3 r t . (3) 

Using reduced coordinates corresponds to switching to a scaled coordinate system where the density is unity. We term 
a microscopic configuration physically relevant with respect to a given thermodynamic state point if the configuration's 
contribution to the partition function at that state point is not a priori negligible. For instance, no configurations 
where all particles occupy the left half of the system's volume are physically relevant for ordinary liquid states. 

State points (1) and (2) with temperatures T\ and T 2 and densities p\ and pi are isomorphic if they obey the 
following: Whenever two of their physically relevant microscopic configurations (r^\ ... , r^') and (r^\ ... ,r^) have 
identical reduced coordinates (i.e., f ^ = f^ 2 '), they have proportional configurational NVT Boltzmann factors: 

e -U(v { ^\ ... ,r< v 1 »)/fe B T 1 = e -U(r[ 2 \ ... ,rff)/k B T 2 ^ 

It is understood that the constant C12 depends only on the state points (1) and (2), not on the microscopic configu- 
rations. Isomorphic curves in the state diagram are defined as curves on which any two state point are isomorphic. 
An IPL liquid with interactions scaling with distance cx r~ n trivially obeys Eq. Q with C\<x = 1 for states with 

p^ 3 jT\ — p 7 ^ 3 /T2- No other systems obey Eq. Q rigorously, but we show in Sec. Ill from simulations of a strongly 
correlating liquid that the existence of isomorphs is a good approximation. Although only IPL liquids have exact 
isomorphs, we shall say that a liquid "has isomorphs" if isomorphic curves exist to a good approximation in the 
liquid's phase diagram. 

The isomorph definition does not refer to WU correlations, but only strongly correlating liquids have isomorphs. 
The Appendix proves in detail that all strongly correlating liquids have isomorphs and vice versa. It is illuminating 
here to briefly sketch why the existence of isomorphs for a liquid implies that it must be strongly correlating: Consider 
a liquid with two isomorphic state points that arc infinitcsimally close to each other. If 5 represents the variation 
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between two infinitesimally close microscopic configurations of the state points with same reduced coordinates, taking 
the logarithm of Eq. Q leads to 5{U/T) — Const, where the constant is infinitesimal. The differentiation leads to a 
relation of the form SU = (da)U + db where da and db are infinitesimal. Since SU is the potential energy difference 
between two microscopic configurations differing by Sri oc r^, one has SU oc W (Appendix A). Altogether we get 
W = AU + B for some constants A and B. This implies 100% correlation of the WU fluctuations - recall, however, 
that the existence of isomorphs is itself an approximation. Thus liquids having isomorphs to a good approximation 
must be strongly correlating. 

A number of properties characterize isomorphic curves in the state diagram of a strongly correlating liquid. Most 
isomorph properties come in the form of isomorph invariants. These are consequences of the fact that for any two 
isomorphic state points there is a one-to-one correspondence between the state points' physically relevant microscopic 
configurations, such that corresponding configurations have the same NVT canonical probabilities. This one-to-one 
correspondence motivates the name isomorph ("same form"), which is fundamental throughout mathematics. Here 
two objects are termed isomorphic if they are structurally equivalent, i.e., if a structure-preserving bijective mapping 
between them exists. In physics and chemistry isomorphic crystals by definition have symmetry groups that are 
mathematically isomorphic - such crystals have the same structure but different constituents. 



B. Isomorph properties 

Before detailing the consequences of the isomorph definition Eq. p} we briefl y recall well-known facts of the 
statistical mechanics of classical liquids in the NVT canonical ensemble ! 8 ! 44 ! 45 * 4 * ^ 4 ^ The Helmholtz free energy F is 
the sum of an ideal gas term and an "excess" free energy term reflecting the molecular interactions, F — Fa + F cx . 
The first term is the free energy of an ideal gas at the same density and temperature, = — NkgTln^A 3 ) where 
A = h/y/2nmkBT is the thermal de Broglie wavelength. The excess free energy is given by 



-F BX /k B T 



dri drN_ £/( ri ,..., rjv )/fc B T 

v V 



(5) 



The integral involves only the configurational degrees of freedom. This separation of the configurational degrees of 
freedom from the momenta is assumed throughout this paper. Thus when we refer to the canonical ensemble and 
canonical probabilities, only the configurational part of phase space is implied. Likewise, the microcanonical ensemble 
below refers to the uniform probability distribution on the constant potential energy surface in configuration space. 
The configuration space probability distribution normalized to the above dimensionless integral is given by 



P(ri,...,rjv) = e-^'-wW^. (6) 

The excess entropy S cx is defined by S cx — —dF cx /dT. Since the entropy in the canonical ensemble is generally given 
by — fcs(lnP), we have 

S cx = -k B J ^y-...^Y-P(ri,...,r N )lnP(ri,...,r N ) . (7) 

The inequality — PlnP < 1 — P implies that S cx is always negative, a fact that is physically obvious since any liquid 
is more ordered than an ideal gas at the same volume and temperature. 

Most isomorph invariants are consequences of two fundamental isomorph properties, Eqs. Q and ^ below. 
Equation Q implies that the normalized reduced-coordinate probability distribution is invariant along an isomorph: 



P(f i, f at) is an isomorph invariant . (8) 

The notation P is introduced to distinguish from the P of Eq. P is normalized via J P(?i, jv)dri...dr.iy = 1 
which implies P(fi, ...,fjv) = N~ P(r±, ...,rjv). The second fundamental isomorph property is that an isomorph I is 
characterized by two functions, fi(r±, ...,fjy) and g(Q) where Q denotes the state point, such that for any physically 
relevant microscopic configuration of state point Q 



U(r 1 ,...,r N ;p) = k B T //(f i, f N ) + g(Q) 



(9) 
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This follows from the isomorph definition Eq. Q. In this formulation the potential energy function is formally 
regarded as density dependent, a dependence that is to reflect the fact that Eq. ^ only applies for physically 
relevant microscopic configurations, i.e., configurations that fill out the volume. 

Before deriving the isomorph properties we note two multidimensional geometric isomorph characterizations. Recall 
that the potential energy landscape is the graph of the potential energy function, i.e., a subset of R 3N+1 . The first 
geometric property is that isomorphic state points have potential energy landscapes which, when restricted to the 
physically relevant states, are identical except for a vertical displacement and scaling by the inverse temperature, and 
a horizontal scaling to unit density. This follows from the isomorph definition Eq. Q, and this is what Eq. ^ 
expresses. The second geometric isomorph characterization relates to the hypersurface f2 in R 3N where the potential 
energy equals the average potential energy of the state point in question: In reduced coordinates this "constant 
potential energy hypersurface", denoted by f2, is invariant along an isomorph. This follows from the first geometric 
property in conjunction with the isomorph definition Eq. Q. For details, please consult Appendix A, which gives 
detailed proofs that a liquid is strongly correlating if and only if it has isomorphs, and that this happens if and only 
if the liquid has curves in the state diagram along which fi is invariant. 

1. Thermodynamics 

A number of thermodynamic quantities are invariant along isomorphic curves in the state diagram of a strongly 
correlating liquid. 

• la. The excess entropy S cyl is invariant along an isomorph. Equation ^ implies that S C7C = 
— kg J Pin Pdri...drpf + Const., from which property la follows because of the invariance of P. Property 
la may also be derived by referring to the microcanonical ensemble where the excess entropy is ks times the 
logarithm of the area of £1: Because f2 is an isomorph invariant, so is the excess entropy. 

• lb. The configurational entropy 5 con f is invariant along an isomorph. The term "configurational entropy" - not 
to be confused with S c -^ - is used here in the sense of Adam and Gibbs and s ubseque nt workers, who related S^onf 
to the temperature dependence of viscous liquids' average relaxation timeP^IGillll The configurational entropy 
is ks times the logarithm of the "density of states" of pote ntial energy minima, the so-called inherent states, 
evaluated at the state point's average inherent state energyP^22l Property lb follows from the identity of the 
scaled potential energy landscapes of two isomorphic state points. 

• lc. When a liquid is heated slowly along an isomorph, the measured specific heat equals the ideal-gas specific 
heat. Since the specific heat is dS/dhiT and property la implies dS cyi — along an isomorph, the measured 
specific heat equals the specific heat that would be measured for an ideal gas subjected to the same process. 

• Id. The isochoric specific heat is invariant along an isomorph. If the liquid's (extensive) excess isochoric specific 
heat is denoted by C™, in terms of the variable X = U jkoT Einstein's expression = ((AU) 2 } /kgT 2 becomes 
Cy x = fcs((AX) 2 ). Equation (JoJ) implies that for two isomorphic state points their microscopic configurations 
with same reduced positions obey X\ — X^ + Const., i.e., for the fluctuations AAi = AA2. Since AX depends 
only on the reduced coordinates and since P(?i, ...,f jy) is isomorph invariant, it follows that the excess isochoric 
specific heat is an isomorph invariant. This implies invariance of the full C™, because the contribution to 
from the momentum degrees of freedom is state-point independent. 

2. Structure 

Particle distribution functions arc generally invariant along an isomorph when quoted in reduced coordinates. 

• 2a. Scaled radial distribution function(s) - as well as higher-order equilibrium particle probability distributions 
- are invariant along an isomorph. Property 2a follows from Eq. (j8j; note that it applies for liquids with any 
number of different types of particles. 

• 2b. The multiparticle entropies S2, S 3, ... are invariant along an isomorph. The equilibrium particle distributions 
give rise to n-particle entropie^^l contributing to the total entropy as follows 5 CX = S2 + S3 + .... In terms 
of the radial distribution function for a system of identical particles g(r), the pair-correlation contribution S2 
is given by S2/N = —(pks/2) J dr [g(r) In g(r) + 1 — g(r)]. When this expression, as well as the more involved 
expressions defining S3, etc., are rewritten in terms of reduced coordinates, it becomes clear that property 2b 
is a consequence of Eq. (J8j. We considered here only the case of identical particles, but property 2b holds for 
systems with any number of different particles. 
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Quantity 


Newtonian dynamics B 


rownian dynamics 


Energy unit (Eo) 






Length unit (Zo) 


-1/3 


-1/3 


Time unit (to) 


\/m/kaT 


P 2 / 3 pJ= R T 


Diffusion constant: 






D = L>/(l /t ) 


E> = \P y/m/KBJ- ) L) 


^ = pk B T D 


Viscosity: 






fj = v/(E t /l 3 o) 


^ P V^ mkB T 71 




Heat conductivity: 






ii = K/(fes/io*o) 


^/m/k B T 
K ~ P 2 / 3 fc B ^ 





TABLE I: Reduced units. The energy and length units refer to state-point properties only, the time unit depends also on 
the dynamics. Once these three fundamental units have been defined, transport and other properties have uniquely defined 
dimensionless versions, denoted by a tilde. The table gives three examples of such properties. 



3. Equilibrium dynamics 

As for the static isomorph invariants, dynamic invariants also derive from the fact that for all physically relevant 
microscopic configurations the potential energy landscapes of isomorphic state points are identical - except for additive 
constants and overall scalings that do not affect the reduced dynamics. For reference below, Table [T] summarizes the 
definition of the basic units and gives examples of some of the corresponding reduced quantities. 

• 3a. Both NVE and NVT Newtonian dynamics are isomorph invariant when described in reduced units. Consider 
first standard energy-conserving Newtonian dynamics, the NVE ensemble. If the mass of particle i is m«, 
Newton's second law is mjfj = Fj where Fj = —V Ti U is the force on the i'th particle. We rewrite this in 
terms of reduced units as follows. If the average particle mass is m, the reduced mass of the i'th particle is 
defined by fhi — rrii/m. The reduced potential energy U is defined by U = U/ksT and the reduced force 
Fj by Fj = —VfJJ. We also define a reduced time, t = t/t a , where to = p -1 / 3 '^[mfkffl '. In terms of these 
reduced variables Newton's second law becomes rhjf j = Fj. The isomorph invariance now follows from Eq. ([9]), 
which implies that the reduced force is the same function of the reduced particle positions for all state points 
on a given isomo rph. Proceeding to Newtonian dynamics in the NVT ensemble realized via the Nose-Hoover 
thermostat J^^l recall that the forces here have the additional term —C^rriiVi with a "friction constant" £ obeying 
£ = (K/K — 1)/tq, where K is the kinetic energy, K its average, and r the thermostat time constant. These 
equations become isomorph invariant when rewritten in terms of the same reduced units as the NVE Newtonian 
equations, if the thermostat time constant is adjusted to be the same in reduced units; otherwise, only the 
long-time NVT dynamics is isomorph invariant. 

• 3b. Brownian dynamics is isomorph invariant when described in reduced units. The Brownian (Langevin) equa- 
tions of motion are first-order stochastic equations. These equations obey detailed balance, ensuring consistency 
with the canonical ensemble. The Brownian equation of motion is rj = /iFj + £(i) where /x is the "mobility" (ve- 
locity/force) and £(t) is a Gaussian white-noise term characterized by (€m(t)£n(t')) = 2fjbkBT8 mn 8[t—t') (m, n — 
1, 2, 3). The path-probability functional is giverPS by P cx exp[— l/(4/x/cs2 n ) Si /_ — ^i) 2 di\. We rewrite 
this in terms of reduced variables with U — U /ksT as above, but a reduced time that is now defined via t = t/to 
where to — p~ 2 / 3 / [iksT . This leads to P cx exp[— l/4^] i (fj — Fj) 2 dt]. Since the reduced force is isomorph 
invariant for microscopic configurations with same reduced coordinates, it follows that reduced-time Brownian 
dynamics is isomorph invariant. 

• 3c. Normalized time- autocorrelation functions - as well as normalized higher-order time correlation functions 
- are invariant along an isomorph when quoted in reduced units. Consider the time-autocorrelation function or 
higher-order time-correlation functions of some variable A referring to constant volume dynamics. Properties 3a 
and 3b imply that for both Newtonian and Brownian dynamics time-autocorrelation functions of A are invariant 
as functions of the reduced time if they are normalized by dividing by (A 2 ). Normalization of any higher-order 
time-correlation function similarly makes it isomorph invariant as function of the reduced times. 
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• 3d. Average relaxation times are isomorph invariant when quoted in reduced units. For any variable A with zero 
mean a generic definition of its average relaxation time is ta = (A(0)A(t))dt/ (A 2 ) . In reduced units this 
becomes ta = f™(A(0)A(t))di/(A 2 ). By property 3c this expression is isomorph invariant. 

• 3e. Reduced transport coefficients like the diffusion constant, the viscosity, etc., are invariant along an iso- 
morph. By the fluctuation-dissipation (FD) theorem the diffusion constant is given by D = L°° (v x (0)v x (t))dt 
where v x is the x component of a particle's velocity. The reduced diffusion constant D is defined by 
D = (p^ 3 y/m/k B T)D for Newtonian and D = D/(pk B T) for B rownian dynamics (Table Q). In both cases 
one has D = (v x (0)v x (t))dt, implying that D is an isomorph invariant because both dynamics are isomorph 
invariant. Similarly, if rj is the viscosity, the reduced viscosity is defined by 77 = (p -2 ^ 3 / \/rnk B T)r] for New- 
tonian dynamics and r) = (/zp -1 / 3 )?? for Brownian dynamics. When rewritten as a reduced-time integral over 
the reduced shear-stress autocorrelation function, the required result follows. Similar results apply for the heat 
conductivity and other DC transport coefficients. 

• 3f. Goo/Tp is invariant along an isomorph where Goo is the instantaneous shear modulus. If S xy = X^ 3 -*-^*,!/ 
where Xi is the x-coordinate of the i'th particle and F i y is the y-component of the force acting on it, the FD 
expression for the instantaneous shear modulus^ is Goo = pk B T + (S xy ) /Vk B T . In terms of reduced variables 
one has XiF^ y = —k B TxidU / 'djji, leading to Goo/pk B T = 1 + {(^ Ji x i dU / dyA 2 ) / N . The required isomorph 
invariance now follows from Eq. ^ . 

4. Aging 

Not all isomorph properties come in the form of invariants. This section discusses a different type of consequence 
of the existence of isomorphs for strongly correlating liquids. It is assumed that the externally controlled variables 
are volume and temperature. 

• 4a. A jump between two isomorphic state points starting from equilibrium takes the system instantaneously to 
equilibrium. This is because the normalized Boltzmann probability factors are identical for the two systems. 
Thus isomorphs are predicted to be a kind of "wormholes" in the state diagram along which one can jump 
instantaneously from equilibrium to equilibrium, even when the states are characterized by long relaxation 
times. 

• 4b. Jumps from any two isomorphic equilibrium state points to a third state point lead to the same aging behavior 
for all physical quantities. Note first that a jump between two arbitrary state points starting from equilibrium, 
1^3, has the same relaxation pattern as the 2 — > 3 jump, where state point (2) is isomorphic with state point 
(1) and has the same density as state point (3): Suppose that instead of the 1 — * 3 jump we first impose the 
isomorphic jump to state point (2) and then, immediately thereafter, jump to state point (3). On the one hand, 
the system will never "register" it spent a tiny amount of time at state point (2). On the other hand, the 1 — ► 2 
jump took the system instantaneously to equilibrium at state point (2) (property 4a). Consequently, the 1 — > 3 
and 2 — > 3 jumps must have the same relaxation towards equilibrium for all physical quantities. It now follows 
that if state point (1) is replaced by an isomorphic state point (1'), the 1 — ► 3 and 1' — > 3 relaxations towards 
equilibrium are identical. 

A concise way of summarizing the aging properties of strongly correlating liquids is that isomorphic state points 
are equivalent during any aging scheme. 

C. IPL invariants and general isomorph invariants 

As mentioned, an IPL liquid has exact isomorphs. These are the curves in the state diagram given by the equation 
p n / 3 /T = Const. All above isomorph properties apply exactly to IPL liquids (single- or multicomponent, as long as 
the exponent is the same for all particles) . Only some of the IPL invariants along the curves given by p"/ 3 /T = Const., 
however, give rise to general isomorph invariants. Examples of IPL invariants that do not generalize are: F cx /T, U/T, 
W/T, the excess pressure coefficient [3y — (d(W/V) /dT) v , and K T /Tp where K T is the isothermal bulk modulus. 
One way of determining which IPL thermodynamic properties give rise to general isomorph invariants is to use the 
approximate equation for the excess free energy derived in Paper IIlP Here it was shown that the hidden scale 
invariance of a strongly correlating liquid implies that one may write F ex (V,T) — f(V) + NksT^p 1 /T) to a good 
approximation. For an IPL liquid /(V) = which, however, does not apply generally. The existence of the f(V) 
term implies that IPL properties involving volume derivatives do not give rise to general isomorph invariants. 
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D. Identifying isomorphs 

How to identify the isomorphic curves in a strongly correlating liquid's state diagram? As we just saw, for an 
IPL liquid the answer is simple: All static and dynamic IPL invariants refer to the curves given by p 1 /T = Const, 
where 7 = n/3. If the virial and potential energy equilibrium fluctuations from their average values obey AW — 
jAU, the number 7 may be expressed in terms of equilibrium fluctuation averages in three simple ways: 7 = 
(AWAU)/((AU) 2 ) = {{AW) 2 ) / {{AU) 2 ) = {{AW) 2 )/{AWAU). This applies only for a 100% correlating liquid, 
however, i.e. for an IPL liquid. For a general strongly correlating liquid there are three different corresponding 
gammas, 



(AWACO ((AWT) 73 = ((AVtQ 2 ) 

71 {{AU) 2 ) ' 72 V ((AC/) 2 ) ' 73 {AWAU) ' 1 ' 

It follows from the definition of the correlation coefficient R (Eq. ^) that 

71 = i? 72 = i? 2 73 . (11) 
Thus 72 is the geometric mean of 71 and 73 , and 



7i < 72 < 73 • (12) 

Although for any strongly correlating liquid the three gammas are quite similar, the question is which gamma to 
use to identify the isomorphs? The answer is that there is no unique gamma. To see this, note that in complete 
generality any quantity q of any liquid defines a state-point dependent exponent j q {Q) with the following property: 
An infinitesimal change away from state point Q conserves q whenever the quantity p^(Q) /T is kept constant (thus 
lq{Q) = {d\nT / din p) q {Q)). Strongly correlating liquids are characterized by the particular property that these 7's 
are all very similar for the isomorph invariants. 

In the simulations reported in the next section we used the excess entropy's "density scaling exponent" 7 derived 
as follows. Along a configurational adiabatic curve = dS cx = {dS ex /dV)TdV + {dS cx /dT) v dT. The volume- 
temperature Maxwell relation for the configurational degrees of freedom implies {dS cx /dV)r — [d{W /V) / dT)v ■ 
Thus a configurational adiabat is characterized by {din p){dW/dT) v = {dlnT)T{dS cx /dT) v = {d In T){dU/dT) v , 
i.e., 



dlnT\ (%£) v _ {AWAU) 



This is 71 of Eq. (13 1; the last equality is a standard thermodynamic fluctuation identity derived, e.g., in Appendix 
B of Paper I. Since it is convenient to think of the density scaling exponent as a generic quantity, we will not refer to 
it as 71, but simply as 7. Note that 7 = {dW/dU)v, which implies that 7 is the slope of the lines of average virial 
versus average potential energy at constant density in the UW plot (compare Fig. 4 in Paper I). 

We propose to use the excess entropy's density scaling exponent because, on the one hand, it rigorously reproduces 
property la and, on the other hand, it may be calculated from equilibrium fluctuations at any given state point. It 
may be shown that for a given infinitesimal density change the optimal temperature change leading to a mean-square 
best fit of the logarithm of Eq. Q is obtained by using this 7 (Appendix B). Finally, note that this 7 is also the 
one suggested by taking the ratio of Eqs. (28) and (27) in Paper III: 7 = ffi/c^ where (3^ = {d{W/V)/dT) v is the 
excess pressure coefficient and c^? is the excess isochoric specific heat per unit volume. 

It is now clear how to step- by-step map out an isomorph which is realized as a configurational adiabat: At any 



given state point one calculates 7 from the equilibrium fluctuations using of Eq. ( 13 1. For a slight density change one 
then determines the temperature that keeps p 7 /T unchanged from the initial state point. This identifies a new state 
point that is isomorphic with the initial one. At the new state point a new value of 7 is calculated, etc. We used this 
method in the simulations reported below, but found almost as good agreement with predictions using 72 (73 was not 
tested). 

In the below simulations, as well as in previous simulations, we generally found only weak state-point dependence 
of the density scaling exponent 7. As shown in Appendix B, for a strongly correlating liquid any variation of 7 with 
state point can only come via a density dependence: 7 = "f{p). This result is derived in the "isomorph approximation" 
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FIG. 1: Direct check of the isomorph condition Eq. Q for the Kob- Andersen binary Lennard-Jones (KABLJ) liquid with 8000 particles. 
The molecular dynamics algorithm generates a time sequence of equilibrium microscopic configurations at the state point (1) given by 
pi = 1.258 and T\ = 0.628 in standard LJ units. For each configuration the potential energy is plotted on the x-axis, while the y-axis 
gives the potential energy of the same microscopic configuration scaled to density p2 = 1.228. Clearly these two quantities are highly 
correlated, which is the requirement for an isomorph. The slope of 0.884 implies that state point (P21T2) is isomorphic with state point 
[pi, Ti) if T2 = 0.884Ti = 0.555. In comparison, the exponent of Eq. evaluated from the fluctuations at state point (1) is 7 = 5.018, 
and the predicted isomorph temperature at state point (2) calculated from the requirement of keeping p 7 /T constant is T2 = 0.556. The 
non-zero offset of the best fit line (of —0.909) reflects the fact that the constant C12 in Eq. (Ul differs from unity. 



where there is identity of the curves in the phase diagram of constant excess entropy and those of constant isochoric 
specific heat. Our simulations are consistent with this result in the sense that 7 varies significantly less along isochores 
than otherwise throughout the state diagram. The conclusion is that, in the isomorph approximation, any isomorph 
invariant is a function of p 7 ^/T, i.e., it may be written as a function of a variable of the form e(p)/T. This 
is the original recipe for collapsing measurements at different pressures and temperatures, which was proposed by 
Alba-Simionesco and her coworkers in their pioneering paper on density scaling! 2 ^ 



III. COMPUTER SIMULATIONS 



This section presents results from computer simulations investigating some of the predicted isomorph properties. 
The purpose is to document the existence of isomorphs for a typical strongly correlating liquid. Results are reported 
for molecular dynamics simulations of the standard Kob- Andersen 55 56 80:20 binary Lennard-Jones liquid (KABLJ) 
with N = 8000 particles simulated using the Gromacs 

package. 1 ^ 1 ] The KABLJ system is a strongly correlating 
liquid 1 3 which is easily supercooled without crystallizing. 59 Its density scaling exponent 7 varies slightly with state 
point, but at low and moderate pressure and temperature 7 stays between 5 and 6 (at extremely high pressure and 
temperature it converges to the value 4, which a purely repulsive r -12 potential would imply). The simulations 
were performed in the NVT ensemble using the Nose-Hoover thermostat with characteristic time 0.5 in MD units. 
Collections of isomorphic state points were identified as described above. 



A. Direct isomorph check 



It is possible to check directly the proportionality between Boltzmann factors which defines an isomorph. This is 
done as follows. At one state point of the equilibrium KABLJ liquid the simulation generates a time sequence of 
microscopic configurations. We now ask: For a given density change, does a temperature exist at which the new 
density state point is isomorphic to the initial one? If yes, what is this temperature? These questions are answered by 
plotting the potential energy of each microscopic configuration against the potential energy of the same configuration 
scaled to the new density. In Fig. [I] density was decreased by 2.4%, which corresponds to a decrease of the relaxation 
time by more than a factor of four if temperature is not changed. The two potential energies are 99.9% correlated. This 
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FIG. 2: (a) AA particle radial distribution functions for the KABLJ liquid at six isomorphic state points, (b) The same radial 
distribution functions plotted as functions of the reduced distance, showing very good data collapse. The inset focuses on the 
first peak and includes results for a second collection of isomorphic state points, (c) AB and BB radial distribution functions 
in reduced units for the six isomorphic state points of (a). There is good data collapse, but with larger deviations than for the 
AA distribution function. This shows that isomorph properties are not exact. 



shows that a temperature does exist at which the isomorph condition Eq. Q is fulfilled to a very good approximation. 
When the density change goes to zero, the correlation goes to 100% of course, but the high correlation is nevertheless 
noteworthy given the fact that the slope of the stretched oval in the figure is not very close to unity (the slope is 
0.884). The interpretation of the slope is that the two state points are isomorphic if the temperature at the new 
density is 0.884 times the old temperature. - Appendix C details how the correlation coefficient of Fig. [T] relates to 
the standard WU correlation coefficient of Eq. ([2| . 

One might think that Fig. [IJproves that all the isomorph invariants apply to a good approximation for the KABLJ 
liquid, but this is not necessarily the case. Consider for instance the dynamics. This becomes increasingly barrier 
dominated as temperature is lowered,^ and the increasingly unlikely event of a barrier transition determines the 
relaxation time. Even if almost all microscopic configurations of significant weight in the canonical ensemble obey the 
isomorph definition Eq. Q very well, this does not ensure that the barriers scale in the same way. Consequently we 
need more simulations to validate the isomorph concept. 



B. Equilibrium properties: Statics and dynamics 

Denoting the large Lennard- Jones particle as A, Fig. [2]ja) gives the AA radial distribution functions for six 
isomorphic state points of the KABLJ system. The temperature varies by almost a factor of two; nevertheless 
there is good collapse of the curves when these are plotted as functions of the reduced distance (Fig. [2](b)) . The 
inset of (b) zooms in on the peak of the radial distribution function. We added here data for a second collection of 
isomorphic state points, showing that different isomorphs have different structure. That the isomorph invariants are, 
after all, not exact is clear from Fig. [2^c), which shows the scaled AB and BB radial distribution functions of the 
isomorphic state points of (a). 

Turning now to the dynamics, Fig. [3ja) shows the AA self part of the intermediate scattering functions F s (q,t) 
at the reduced wavevector corresponding to the first peak of the static structure factor, where t is time scaled 
by the characteristic time to = ^^(m/fesT) 1 / 2 (Table [i]). The figure shows results for four state points at two 
densities. Two of the four state points are isomorphic. These two state points (black and red curves) have virtually 
identical relaxation behavior, including the short-time "cage-rattling". Figure [3] confirms isomorph properties 3c and 
3d. For the two isomorphic state points of Fig. [3]ja), (b) shows the four-point dynamic susceptibility defined^ by 
X4(t) = N((AF s (q, i)) 2 ) , i.e., the mean-square fluctuation of the self-part of the intermediate scattering functions. 
The quantity Xi(t) measures the degree of dynamic heterogeneity on a given time scale. By the predicted isomorph 
invariance of time-correlation functions as functions of reduced time (property 3c) X4W is predicted to be invariant, 
which is indeed the caseP^Ml 

As shown briefly in Sec. II and in detail in the Appendix A, only strongly correlating liquids have isomorphs. 
A system that is not strongly correlating^ is the SPC water model, where the hydrogen bonds are mimicked using 
Coulomb interactions and the oxygen atoms interact via an LJ potential.^ The fact that this model has near zero WU 
correlation reflects^ water's density maximum. In order to prove that SPC water does not have isomorphs, suppose 
that it did. Then state points with same reduced diffusion constant would have same structure. Figure |4]ja) shows 
the reduced diffusion constant as a function of temperature for two sets of isochoric state points. Interpolation with 
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FIG. 3: (a) The self part of the intermediate scattering functions at the wavevector corresponding to the first peak of the static structure 
factor as function of the reduced time t for four state points of the KABLJ liquid. Two of the state points are isomorphic (red and black 
curves); these two state points have virtually indistinguishable self part of the intermediate scattering functions. 

(b) The four-point dynamic susceptibility X&W f° r the two isomorphic state points of (a). This quantity, which measures the degree of 
dynamic heterogeneity, is predicted to be an isomorph invariant. The full curve is a cubic spline fit to the black data points. 




FIG. 4: (a) Reduced diffusion constants for two sets of isochoric state points of the SPC water model, which is not strongly 
correlating. The data were fitted by a polynomial in order to identify the temperature where the two densities have the same 
reduced diffusion constant (T=252K). 

(b) Oxygen-oxygen radial distribution functions as functions of reduced distance for the two state points of SPC water identified 
in (a) with same temperature and same reduced diffusion constant, but different density. If water had isomorphs, the two radial 
distribution functions would be virtually identical. 



a polynomial was done in order to identify the temperature where the two densities have the same reduced diffusion 
constant (T=252 K). Figure |4jb) shows the radial distribution functions of these two state points. Clearly, SPC water 
does not have isomorphs. 



C. Out-of-equilibrium properties: Aging 

We showed by example that isomorphic state points have the same scaled static and dynamic correlation functions, 
but all properties tested so far were equilibrium properties. What happens when a strongly correlating liquid is taken 
out of equilibrium? To answer to this we simulated temperature/density jumps from equilibrium, i.e., instantaneous 
changes of these two variables to new values. All states involved in the "aging experiments" belong to the two sets 
of isochoric points whose self-intermediate scattering functions and radial distribution functions are plotted in Figs. 
[2j» and [2ft)). 

FigureT5Ta) shows the time evolution of the potential energy when a jump is made from equilibrium, bringing 
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FIG. 5: (a) Results from simulating an instantaneous temperature and density jump applied to the KABLJ liquid, starting 
from equilibrium and jumping to a state which is isomorphic with the initial state (LJ units). Except for the transient no 
relaxation is associated with the jump, showing that the system is instantaneously in equilibrium. The relaxation time of 
both state points is about 500 in reduced units (the two state points have the same reduced relaxation time because they are 
isomorphic). 

(b) Potential energy relaxation towards equilibrium for the KABLJ liquid for different jumps: 1 — ► 2, 3 — ► 2, and 4-^2 (see 
the inset). The 1 — > 2 jump is the isomorphic jump of (a), the two other jumps are not between isomorphic state points. - 
The results shown in (a) and (b) are averaged over ten independent simulations. 
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FIG. 6: Potential energy relaxation towards equilibrium for the KABLJ liquid comparing the jumps 1 — > 3 and 2 — » 3 (see the 
inset). The two relaxations are predicted to be identical because state points (1) and (2) are isomorphic. The results shown 
are averaged over ten independent simulations. 

the KABLJ liquid to a new state point which is isomorphic to the initial state. For both state points the reduced 
relaxation time is around 500. The jump was performed as follows: We instantaneously increased the box volume 
without changing particle positions (the initial overshoot is due to this) and simultaneously changed the thermostat 
temperature to the final temperature. There are no signs of slow relaxation after the jump. Thus the system is 
immediately in equilibrium, as predicted for jumps between isomorphic state points (property 4a). 

The horizontal (red) line in Fig. |5jb) shows the data of Fig. (5fa), now on a logarithmic scale with time shifted such 
that the jump occurs at t — 0. Here it is even more clear that jumps between two isomorphic state points preserve 
equilibrium. In contrast, the 3 — > 2 and 4 — > 2 jumps both age slowly to equilibrium, where all three initial states 
(1), (3), and (4) have the same density. Clearly, instantaneous equilibration is a feature of jumps between isomorphic 
points only. 
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The preservation of equilibrium for jumps between isomorphic state points has important consequences for general 
jumps (Fig. [BJ. Consider first the inset of Fig. [6j Suppose we start in equilibrium at state point (1) and change 
temperature and density to T 3 and p%. According to property 4b, if state point (2) is isomorphic to state point (1) 
the observed relaxation behavior for the 2 — > 3 jump should be the same as for the 1 — > 3 jump. The simulations 
confirm this prediction. 

We can now also understand why the 3 — > 2 jump in Fig. [5^b) approaches the equilibrium potential energy from 
below, although state point (3) has an average potential energy which is slightly higher than that of state point (2) 
(compare the inset of Fig. [5](b)) . This is because, whenever a jump is performed, the system makes an instantaneous 
isomorphic jump to the new density to relax from here to the final state point by moving on the isochore. Thus in 
Fig. |5^b), the 3 — > 2 jump is first an isomorphic jump to the state with correct density, a state that has lower average 
potential energy than state (2). Only thereafter the system relaxes towards state (2). 

In most experiments pressure, not density, is controlled. The only difference is that an instantaneous isomorphic 
jump goes to a new pressure, instead of to a new density, and that the subsequent relaxation follows an isobar instead 
of an isochore. Thus relabeling the dashed lines in the insets of Figs. [BJb) and [6] to isobars would qualitatively 
illustrate an experimental temperature and pressure jump. 

IV. RELATING ISOMORPHS TO SOME TOPICS OF CURRENT LIQUID STATE THEORY AND 

EXPERIMENT 

This section discusses connections between the isomorph concept and current liquid state theory and experiment. 
No subject is treated in depth; the purpose is merely to show by example that isomorphs fit nicely into a number of 
previous findings. In most cases we demonstrate how these are consistent with the isomorph concept, but in a few 
cases a previous finding is shown to be a consequence of the existence of isomorphs for strongly correlating liquids. 

A. Phenomenological melting rules 

Melting is an old subject of condensed-matter physics. There has long been a good understanding of the statistical 
mechanics of melting via the density functional theory of Ramakrishnan and Yussouff! 5 Supplementing this are a 
number of phenomenological melting rules. Perhaps the most famous is the Lindemann criterion, according to which 
melting takes place when the crystal's vibrational root mean-square displacement is about 10% of the nearest-neighbor 
distanced 

For a strongly correlating liquid the melting curve in the phase diagram must be parallel to nearby liquid and 
crystalline isomorphs (the isomorph concept applies to the crystalline phase as well as to the liquid). This is because, 
if an isomorph were to cross the melting curve, the Boltzmann factors favoring crystalline order would dominate 
one part of that isomorph and be negligible on another part; this contradicts the structure invariance. For strongly 
correlating liquids the invariance of the crystalline excess entropy along the melting line implies, as is easy to show, 
pressure invariance of the Lindemann melting criterion. Moreover, since both crystalline and liquid excess entropy 
are invariant along the melting curve, the melting entropy must be pressure independent. 

There are other consequences of melting curves being isomorphs for strongly correlating liquids. Thus along the 
melting curve, slightly to the liquid side, a number of properties are predicted to be invariant: the reduced viscosity, the 
reduced surface tension, the reduced diffusion constant, the reduced heat conductivit y - in fact all iso morph invariants. 
There are several theoretical and experimental works pointing in these directions! 67 * 68 * 69 * 70 * 71 * ^ For instance, the 
reduced-unit static structure factor of liquid iron measured by x-ray scattering is invariant along the melting curve 
for pressures up to 58 GPa. 73 Finally, note that the Hansen- Verlet criterion,^ according to which melting takes place 
when the liquid structure factor peak is 2.85, due to the invariance of the structure factor in reduced coordinates is 
also consistent with the melting curve being an isomorph for strongly correlating liquids. 

B. Rosenfeld's excess entropy scaling 

Rosenfeld presented in 1977 an interesting observation.* 7 ^* 7 ^* For a large class of model systems the reduced transport 
coefficients (table I) appear to be functions of the excess entropy only. He justified what is now known as "excess- 
entropy scaling" by reference to variational hard-sphere thermodynamic perturbation theory.* 7 ^* He emphasized that 
excess-entropy scaling is a semi-quantitative model rather than a theory and referred to excess-entropy scaling as a 
principle of corresponding states. For any given strongly correlating liquid excess-entropy scaling follows from the 
isomorph properties derived in Sec. II, although this does not explain why the functional excess entropy dependence 
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appears to be quasiuniversal.^^] Recent simulations of the Gaussian core model by Truskett and co-workers^ provide 
an example where Rosenfeld's excess-entropy scaling fails. This is not surprising, since this model due to the soft core 
is not expected to be strongly correlating. 

A scaling procedure that is similar in spirit to Rosenfeld's excess entropy scaling was discussed by Dzugutov in 
1996.' 8 He showed from simulations that the reduced diffusion constant follows D oc cxp^So/iVA's) for a number of 
systems, where S2 is the two-particle entropy. This relation is also consistent with the isomorph concept, because 
both D and S2 are both isomorph invariants (properties 3e and 2b). 

From the isomorph viewpoint there is no reason to expect Rosenfeld's or Dzugutov's equations to hold for liquids 
that are not strongly correlating. While genuine Rosenfeld scaling seems to fail for liquids that are not strongly 
correlating, it appears that related excess entropy scaling procedures hold for many such liquids (see, e.g., Refs. l79l 
andlStj]). Why this is, is an important question for future research. 

C. Young and Andersen's approximate scaling principle 

In 2003 Young and Andersen conjectured an approximate scaling principle according to which two liquids have 
similar dynamical properties if they at two same-density state points hav e sim ilar static pair-correlation functions 
- even if their potentials are quite different and the temperatures differJ2S£2] This was confirmed by simulations 
comparing the LJ liquid to the same system with a Weeks-Chandler- Andersen cut-off where only the repulsive part 
of the potential is kept. At same density the latter system required lower temperature in order to have almost 
identical pair-correlation functions. When temperature was properly adjusted, however, the two liquids were shown 
to have similar coherent and incoherent intermediate scattering functions, as well as similar velocity and current 
auto-correlation functions. This confirms the proposed approximate scaling principle for the single-component LJ 
liquid (in contrast, Berthier and Tarjus recently showed that the WCA approximation applied to the KABLJ liquid 
does not reproduce this liquid's dynamics^). 

As pointed out by Young and Andersen, if not only the pair- but all higher-order correlation functions are identical 
for two liquids, they must have proportional NVT Boltzmann statistical weights. This suggests generalizing the 
isomorph concept to an equivalence relation between different liquids. One may define two liquids to be isomorphic 
if: 1) They are both strongly correlating, and 2) a pair of state points exists such that Eq. Q is obeyed for all 
physically relevant microscopic configurations with the same reduced coordinates in the form generalized to two 
different potential energy functions, U\ and \J%. Clearly when this is obeyed, the entire system of isomorphic curves 
of one liquid maps onto that of the other liquid. 

D. Two-order parameter maps of Debenedetti and coworkers 

In a series of publications Debenedetti and collaborators studied liquid structure in terms of a translational order 
parameter, t, and a system-dependent orientational order parameter that may be, for instance, the often used quantity 
q 6 ) 84 | 85 | 86 | 87 | 88 | 89 ] -pjjg translational order parameter is defined as t — J^" \g(f) — l\df where g(r) is the reduced- 
coordinate pair distribution function and f c = p x ^r c is a cut-off; Q§ is defined from sixth-order spherical harmonics 
involving the "bonds" to the twelve particles nearest a given particle (averaged over all particles). The general 
picture one finds by plotting t versus Q 6 in so-called ordering maps is that for the LJ liquid and related systems 
there is a striking collapse, showing that "bond-orientational orde r and translational order are not independent for 
simple spherically symmetric systems at equilibrium" pSl For water^lsU anc j gilicaj^ on the other hand, there is no 
such collapse; here the order parameters cover a two-dimensional region of order-parameter space. These results 
may be explained by reference to isomorphs: Since both order parameters are isomorph invariants, they cannot vary 
independently for a strongly correlating liquid. Water and silica are not strongly correlating (Paper I), so there is no 
reason to expect their two order parameters to follow each other. 

E. Viscous liquid dynamics 

The dynamics of viscous liquids is a long-time focus of our research at the DNRF "Glass and Time" center, and 
the class of str ongly c orrelating liquids was identified in our efforts to understand the properties of single-parameter 
viscous liquids p3 22 l 23 l As shown below, the isomorph concept throws new light on facts and puzzles of the physics of 
viscous liquids. 

1) Cause of the non-Arrhenius relaxation time: The "isomorph filter". A major puzzle concerning viscous liquids 
approaching the glass transition is the origin of the non-Arrhenius temperature dependence of their average relaxation 
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time T l 10 l 12 l 13 l 14 l 15 l 16 l 17 l 18 l 19 l 2 l 21 l 9 ° l 91 l 9: H i n many cases the relaxation-time increase is truly dramatic, with a factor 
of ten or more slowing down if temperature is decreased by just 1%. There are several models for this. Although 
it is not obvious that any simple universal model exists, most workers in the field assumes this to be the case. Any 
universal model must apply to strongly correlating liquids, of course. This implies a criterion which we shall refer to 
as the isomorph filter. The idea is the following. Since the reduced average relaxation time is an isomorph invariant, 
this quantity can only be controlled by a quantity that is also an isomorph invariant. All theories relate lnr to 
some quantity. For any viscous liquid lnr = lnf, because the temperature variation of the average relaxation time 
dominates completely over the reduced time unit's VT temperature dependence. Thus the average relaxation time 
must be controlled by an isomorph inv ariant . Consider now some examples. 

According to the Adam-Gibbs modeP^H the relaxation time r varies as lnr cx l/TS con i(p,T), where S con {(p,T) 
is the configurational entropy discussed in connection with isomorph property lb. Adam and Gibbs wrote the 
proportionality constant as the product of a critical configurational entropy s* and a term Ap which is "largely the 
potential energy hindering the cooperative rearrangement per monomer segment" , and they argued that the state- 
point dependence of these terms can be neglected. This is also what is usually done in comparing to experiments. In 
this case, however, since f and S con {(p,T) are both invariant along an isomorph, the model is inconsistent with the 
isomorph invariance due the extra factor T. Thus in this version the model cannot apply to strongly correlating liquids. 
Alternatively, the proportionality constant must be allowed to vary with state point. Indeed, simulations of Sastry^l 
of the KABLJ liquid show good agreement with the Adam-Gibbs model if the constant is density dependent, a density 
dependence which, in order not to violate isomorph invariance, must be - and to a reasonable approximation is - 
given by a factor proportional to p 1 . Although the original Adam-Gibbs model cannot apply for strongly correlating 
liquids (van der Waals liquids and metallic liquids), it should be noted that for very frag ile liquids the temperature 
dependence of S con f(p,T) dominates over the T factor in the Adam-Gibbs expression pS for these liquids the model 
may work fairly well. Another possibility is that S con {(p,T) alone controls the relaxation time. Interestingly, this was 
suggested by Bestul and Chang in a paper — preceding Adam-Gibbs by one year, which showed data consistent with 
a universal value of the excess entropy at the glass transition for several liquids. Much more recently, Truskett and 
collaborators^ suggested that the excess entropy S ex controls the relaxation time in the spirit of Rosenfeld scalingpilZSI 
this is also consistent with isomorph invariance. As a final example of a theory where entropy controls the relaxation 
time, consider the random first-order transition (RFOT) of Wolynes and co-workers (reviewed in Ref. I57)l . According 
this theory and its generalization by Bouchaud and BiroliJ^for some exponents x and y the relaxation time is given 
lnr cx (Y/T) x S c J nl (p,T) where Y is a surface tension at the molecular scale. This expression is isomorph invariant 
only if Y cx T. This is, however, precisely what was predicted in 2000 on different grounds by Xia and WolynespS so 
the RFOT passes the isomorph filter. 

The free- volume model of Cohen and Grest predicts that lnr cx l/vf(p,T) where Vf(p 7 T) is the free volume. 
As it stands, this model does not survive the isomorph filter because the proportionality constant is state-point 
independent and, although the free volume definition is not obvious, Vf(p,T) is hardly an isomorph invariant. If one 
imagines Vf to be a geometrically determined quantity measured in units of l/p, however, the model is invariant 
along an isomorph. Thus, w hile Vf(p,T) cannot determine the relaxation time, the quantity Vf(p,T)p can possibly. 
The shoving mode l 21 ' 100 ' 1 -^ predicts that lnr cx G oc (p, T)V c /k B T, where V c is a characteristic volume which in 
experiments is of order the molecular volume. This model, which is one of several related "elastic" models,^ is 
consistent with the existence of isomorphs if V c is geometrically determined: In this case V c cx l/p and, because 
Goo(p,T)/Tp is an isomorph invariant (property 3f), this implies that G oa (p,T)V c /kBT is isomorph invariant. The 
vibrational mean-square displacement v ersio n of the elastic models, lnr cx a 2 /(x 2 ), is not isomorph invariant if a 
is as sum ed to be independent of density,^ 2 but it does survive the isomorph filter if the reasonable assumption is 
made 103 that a cx p ' . As a final example, we note that the entropic barrier hopping theory of Schweizer and 
coworker d 104 * 105 * in a certain limit predicts that lnr cx d 2 (F 2 ) / \ksT) 2 where d is an effective hard sphere diameter 
and (F 2 ) is the short-time-averaged single-particle mean-square force. This is expression is not isomorph invariant, 
but the modest modification of it obtained by replacing d by p -1 / 3 leads to an isomorph invariant expression, as is 
easily shown by use ofEq. ^. 

2) Isochronal superposition. Building on earlier works by Roland et aZ.p^in 2005 Ngai, Casalini, Capaccioli, Paluch, 
and Roland published a paper entitled "Do theories of the glass transiti on, in which the structural relaxation time 
does not define the dispersion of the structural relaxation, need revision?"E23 The authors showed here that for many 
viscous liquids and polymeric systems the average relaxation time determines the shape of the dielectric loss peak. 
Thus whether the average relaxation time is increased by lowering temperature or by increasing pressure, the effect 
is the same on the relaxation time distribution as monitored via the dielectric loss. This result, which we for brevity 
refer to as "isochronal superposition", was puzzling at the time (at least to us). Isochronal superposition now appears 
as a consequence of the existence and properties of isomorphs: Any strongly correlating liquid must obey isochronal 
superposition because, according to properties 3c and 3d, both the average relaxation time and the dielectric spectrum 
are isomorph invariants (when given in reduced units, but as mentioned their use makes little difference for viscous 
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liquids). Thus whether temperature is lowered or pressure is increased (or a combination thereof) to bring the liquid 
to state points with the same average relaxation time, all such state points are isomorphic. Consequently, they have 
the same relaxation spectra for quantities probed by linear response experiments. In Ref. 107 hydrogen-bonding 
liquids were quoted explicitly as exceptions to isochronal superposition; this is consistent with the fact that these 
liquids are not strongly correlating.^ 

3) How many parameters are needed for describing a viscous liquid? In the old days of glass science it was actively 
debated whether one or more "order parameters" are required to describe glass structure and the glass transition. 
We write "order parameters" in quotation marks in order to emphasize that this term has a somewhat different 
use in traditional glass science than in the theory of critical phenomena. The original considerations of Prigogine, 
Defay, Davies, and Jones , and others, referred to the glass transition as a second-order phase transition in the 
Ehrenfest sense E5 5 | i09 | xxo | j^. wag recen ^iy shown that strongly correlating liquids are precisely the liquids that to 
a good approximation may be regarded as single-order-parameter liquids a review of the connection was given 
in Ref. [33] Strongly correlating liquids have isomorphs, so it should be possible to link isomorphs directly to the 
singlc-order-parameter scenario. Indeed, by labelling the collection of isomorphic curves in the phase diagram with 
a continuously varying real number, a formal single-order-parameter description throughout the phase diagram is 
arrived at. This suggests a generic way of defining the single-order-parameter scenario that is implicit in the old 
works. The isomorph label is not unique, obviously; it may be chosen as any property that is isomorph invariant, for 
instance C?/. 



F. Some further potential isomorph connections 

1) In 1989 Baranya and Evans^l reported from simulations that the excess entropy of the LJ liquid to a good 
approximation equals the two-particle entropy S2 plus a constant. This is consistent with the existence of isomorphs, 
because on an isomorph both S2 and S'ex are invariant. That is, because the LJ liquid is strongly correlating, if for 
instance S cx were an isomorph invariant but S? were not, the Baranya-Evans finding could not be correct. 

2) Saija, Prestipino, and Giaquinta in 2001 showed by simulations of both the Yukawa and the LJ liquids that the 
scaled radial distribution functions are identical for states where the excess entropy equals the two-particle entropy.E^] 
This is consistent with the fact that these liquids are both strongly correlating, because scaled radial distribution 
function, excess entropy, and two-particle entropy are all isomorph invariants. 

3) The basic idea of mode-coupling theory - that statics determines the dynamics - is consistent with the existence 
of isomorphs: For any strongly correlating liquid one expects that, if two state points have the same reduced pair- 
correlation functions, they are isomorphic. This means that they have the same (reduced) dynamics, so in this sense 
the pair correlation function "determines" the dynamics. From the isomorph perspective mode-coupling theory may 
be expected to work best for strongly correlating liquids. 

4) Recently, Roland briefly reviewed characteristic times and their inv ariance to thermodynamic conditions from 
a general point of view, including also systems that are not liquidsP-^ He showed that quite different transitions 
in systems with slow relaxations - onset of activated dynamics, dynamic crossover in viscous liquids and polymers, 
order-disorder transitions in liquid crystals, vitrification - at varying tempe ratu re and pressure all take place at state 
points with the same value of the relaxation time. The conclusion is thatP»^ "the control parameter driving these 
transitions has the same functional dependence on T, p, and V as the relaxation time." This follows from the existence 
of isomorphs, although we can only explain this observation for strongly correlating liquids and solids. 

5) The concept of hidden scale invariance may have implications beyond liquid state theory. Thus recently Procaccia 
and collaborators studied by simulation plastic flow of amorphous solids in the athermal limit P23 For two-dimensional 
solids composed of multi-disperse particles with the interaction length taken from a Gaussian distribution they showed 
that stress-strain curves at different densities collapse to a master curve. This happens when stress is scaled by p" ', 
where v is 5.87 when the repulsive part of the potential can be fitted by the IPL term r~ 10 . The number 5.87 is not 
far from the 5 predicted by the two-dimensional exponent 7 = n/2 for IPL potentials. 

V. CONCLUDING REMARKS 

This paper introduces the concept of isomorphic curves in the state diagram of a strongly correlating liquid. The 
existence of isomorphs reflects the liquid's hidden scale invariance, and therefore the class of strongly correlating 
liquids is identical to the class of liquids with isomorphs (Appendix A) . Isomorphs may be labelled by any of these 
invariants, for instance Cy~. The isomorph concept, in fact, may be justified by starting from the concept of an "order 
parameter" labelling curves in the state diagram of same physical properties: Suppose a liquid is described by a single 
"order parameter" in the sense that many of its properties are determined by the value of the "order parameter" . In 
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this case one may draw curves in the state diagram along which these properties are all invariant. If one asks how it 
is that several properties could possibly correlate in this way, the simplest answer is that for any two state points on 
a constant "order parameter" curve, the canonical probabilities are proportional for microscopic configurations that 
somehow correspond to each other - and the simplest possibility is that the "corresponding" configurations are those 
that trivially scale into one another. This is nothing but the isomorph definition. Thus the isomorph concept may 
be arrived at via some extrapolation by postulating that single-order-parameter liquids exist and inquiring into their 
properties. 

All isomorph invariants apply to IPL liquids. But as emphasized repeatedly, the converse is not true: Not all 
properties that are invariant for an IPL liquid along states obeying p"/ 3 /T = Const, generalize to become isomorph 
invariants. The exceptions from approximate IPL behavior derive from the fact that the constant C12 of Eq. Q 
generally differs from unity. Only IPL invariants that are independent of the identity C'12 — L which characterizes 
IPL liquids, generalize to the class of strongly correlating liquids. 

The present and previous papers on strongly correlating liquids show that this class of liquids is simpler than 
liquids in general. This is consistent with the general understanding among liquid-state specialists for many years, 
according to which non-associated liquids are generally simpler than associated ones. The virial / potential energy 
correlation coefficient R provides a quantitative criterion for distinguishing simple liquids from the more complex - 
and admittedly often more spectacular - liquids that are not strongly correlating. 
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APPENDIX A: EQUIVALENCE OF THREE CONDITIONS 

This Appendix proves that the following three conditions are equivalent for any liquid: 

• (a) The liquid is strongly correlating. 

• (b) The liquid has isomorphs. 

• (c) The liquid has curves in the state diagram along which the reduced-coordinate constant potential energy 
hypersurface is invariant. 

Only true IPL liquids obey these conditions rigorously (i.e., with 100% correlation in (a)), and the equivalences 
may be stated more accurately as follows: A liquid is strongly correlating if and only if it has isomorphs to a good 
approximation, which happens if and only if the liquid has curves in the state diagram along which the reduced- 
coordinate constant potential energy hypersurface is almost invariant. It is instructive to prove all six implications, 
although logically (a) (6) (c) =$> (a) would suffice. 

Before proving the equivalences we note a few facts. If S denotes the variation between two infinitcsimally close 
microscopic configurations with same reduced coordinates, we first show that 



SU = (dkip) W. (Al) 

For two infinitcsimally close microscopic configurations with the same reduced coordinates (f,; = p 1 / 3 ! - .;) one has 
= Sr i = (p~ 2 / 3 /3)(dp)ri + p 1 ' s Sri, i.e., Sri = — (l/3)(dlnp)r.i. The change of potential energy between the two 
configurations is given by SU = J2i ' V r4 U, which via the definition of the virial, W = —1/3 ^ i • V ri U, leads to 
Eq. (Al). We note further the following identity 



(dlnp) W- (dlnT)U 
f 



(A2) 
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This follows by differentiation and subsequent use of Eq. (All. Finally, note that the infinitesimal version of the 
isomorph condition Eq. Q is 

S Hpj = Const. (A3) 

At any given state point Q, if the average potential energy is (U)q, we define the constant potential energy hypersurface 
f2 as the subset of R 3N given by 

n = {(r u ...,r N )ER 3N \U(r l7 ...,v N ) = (U) Q }. (A4) 
The corresponding reduced-coordinate constant potential energy hypersurface f2 is given by 

Cl = {(fx, r N ) € R 3N | [/(p" 1 / 3 ?!, p-^ 3 r N ) = (U) Q } . (A5) 

Proof that (a) O- (6): A strongly correlating liquid has (near) proportionality between virial and potential energy 
fluctuations, AW — jAU. Thus at any given state point W = 'yll + C to a good approximation for the physically 
relevant microscopic configurations. If density and temperature are changed infinitesimally such that dlnT = jdlap, 



Eq. (A2) implies that 5(U/T) — Const., which is the isomorph condition Eq. (A3 1. Suppose conversely that a 



liquid has isomorphs, and let A denote the difference between two arbitrary, physically relevant configurations at 



the state point in question. Then Eq. (A3) implies AS(U/T) = where 5 (as usual) refers to changes from one 



configuration to another infinitesimally close by with the same reduced coordinates. Via Eq. (A2| this implies 
(dlnp) AW = (dlnT) AC/, i.e., the liquid is strongly correlating with 7 — dhxT/dlnp. 

Proof that (a) <^> (c): For a strongly correlating liquid AW oc AU for fluctuations between physically relevant 
configurations at any given state point. This implies that the hypersurfaces of constant virial and constant potential 
energy coincide. In particular, W is constant on the state point's constant potential energy hypersurface f2. Accord- 
ingly, if density is changed infinitesimally, the change in potential energy between microscopic configurations with the 



same reduced coordinates is the same for all microscopic configurations on £1 (Eq. (Al l). Thus a new hypersurface 
of constant potential energy is arrived at by slightly scaling £1; by adjusting temperature the new hypersurface is 
where the potential energy equals the average potential energy Finally we note that the new and old hypersurfaces 
by construction have same reduced coordinates, thus the two state points have the same f2. Suppose conversely that 
two infinitesimally close state points have the same il. All points on the two constant potential energy hypersurfaces 
differ by the same potential energy, which via Eq. ( |A1[ ) implies that W must be constant on each hypersurface. In 
other words, W is constant on surfaces where U is. This implies 100% correlation between W and U. For large 
systems the fluctuations are small, relatively, and a first order Taylor expansion of the WU relationship leads to a 
linear relationship and 100% correlation between W and U. 

Proof that (b) O- (c): Suppose a liquid has isomorphs. For a state point Q on an isomorph / one concludes from Eq. 
(|9| that (U)q = ksT{fi) + g(Q), where the (canonical) average (//) by Eq. ^ is independent of Q. Consequently, 
= {(?!,.. .,r N ) e R 3N I k B Tf I (r 1 ,...,r N ) + g(Q) = k B T{f z ) + g(Q)} = {(f 1; ...,f w ) e R™ \ /,(fi, ...,f N ) = (/,)} 
is invariant along the isomorph. Suppose conversely that two state points have the same Cl. For these state points 
microcanonical averages are identical for all quantities that may be expressed as functions of the reduced coordinates. 
By the equivalence of the microcanonical and the canonical ensemble, for the two state points in question canonical 
ensemble averages are likewise identical for all quantities that are functions of the reduced coordinates. This can only 
be so if there is identity of the normalized canonical probability factors of any two physically relevant microscopic 
configurations of the two state points with the same reduced coordinates. This is another way of stating that the two 
state points are isomorphic. 

APPENDIX B: PROPERTIES OF THE DENSITY SCALING EXPONENT 

This appendix has two purposes: 1) To derive an optimization property of the density scaling exponent 7 of Eq. 



( 13 1 , and 2) to prove that any possible state-point dependence of 7 comes from a - for strongly correlating liquids 
generally weak - density dependence; more precisely it is shown that 7 = "f{p) in the isomorph approximation. 

As mentioned in the main text there is no unique solution to the problem of finding the "correct" density scaling 
exponent 7, i.e., the exponent identifying isomorphs as the curves along which p 1 /T = Const. First of all, the exponent 
must be expected to vary slightly with state point. But even at a given state point, there is no unique 7 in the sense 
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that all isomorph invariants are mathematically unchanged for infinitesimal steps away from the state point in question 
obeying p 1 /T = Const. This is because, except for IPL liquids, isomorph properties are approximate, so for instance 
the curves of constant excess entropy cannot be expected to be precisely the curves of constant isochoric specific heat 
(or constant reduced relaxation time, etc) - this just applies to a good approximation. 

We argued in Sec. IID that at any given state point there are three obvious gammas, the 71, 72, and 73 of Eq. (10 1. 
We recommend using 71 and used this in Sec. Ill when comparing isomorph predictions to computer simulations. 
This exponent is to be preferred because it - among other things - makes the excess entropy an exact invariant. It 
worked very well with the simulations of the KABLJ liquid, but as noted in the main paper 72 works almost equally 
well. 

The excess entropy gamma 71, henceforth just denoted by 7, has an optimization property coming from answering 
the following question: At any given state point, suppose we change density by an infinitesimal amount. How much 
should temperature be changed to arrive at a new state point which is "as isomorphic as possible" with the original 
state point? To answer this we note that from Eq. (A3 1 one would require the quantity F = S(U /T) to be as constant 
as possible (where the symbol S as previously refers to the difference between two infinitcsimally close microsc opic 
configurations with same reduced coordinates). This is obtained by minimizing ((AT) 2 ). Since according to Eq. (A2| 
one has F oc (dlnp)W — (dlnT)U, for a given density change the quantity to be minimized by varying temperature 

is (((dlnp)AW - (dlnT)AU) 2 ), i.e. 



(din p) 2 ((AW) 2 ) + (dlnT) 2 ((AU) 2 ) - 2(dln p)(d In T)( AW AU) 
Equating to zero the derivative of this expression with respect to <ilnT leads to 



(Bl) 



dlnT _ (AWAU) 
dlnp ~ ((AU) 2 ) 



(B2) 



which is the excess entropy gamma of Eq. ( 13 1 



Suppose instead that one asks the complementary question: For a given an infinitesimal temperature change, what 
is the density change giving a new state point that is as "isomorphic as possible" with the original state point? In 
this case, minimizing Eq. (Bl ) leads to the 73 of Eq. (10 1. Thus if density and temperature were equivalent variables, 



optimization arguments cannot determine which gamma to choose. In the NVT ensemble, however, temperature and 
density are not quite equivalent because, even though they are both externally controlled, the average kinetic energy 
fluctuates while the volume is strictly fixed. 

The second property of the density scaling exponent to be proved is that if the 7 of Eq. (13) varies with state 
point, this variation can come only from a density dependence. More accurately, in the "isomorph approximation" 
where the curves of constant excess entropy and constant excess isochoric specific heat coincide, one has 7 = "f(p). 
Recalling that the excess pressure coefficient (3y is defined by /3y = (d(W/V)/dT)y (Papers I and III) and denoting 
the excess isochoric specific heat per unit volume by dy, the standard fluctuation expressions for these quantities (see, 



e.g., Appendix B of Paper I) implies that Eq. (13 1 may be written 



(B3) 



If there is identity of the curves in the phase diagram of constant excess entropy S ex and those of constant (extensive) 



excess isochoric specific heat C'- 
implies 



(Cy* = Vdy), the definition of the density scaling exponent 7 = (dlnT/dln p)$ c 



dlnT 
dlnV 



(B4) 



From this we get via the mathematical identity (dx / dy) z (dy / dz) x (d z / dx) 



din Cy" 
dlnT 



1 /dlnT 



7 



dlnV 



din Cg? 
dlnT 



1 fdlnC? 



7 



dlnV 



(B5) 



Using Cf? = T(dS cx /dT) v we get via Eq. (B3l 
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dlnC v i \ _ C7 fdlnC^\ _ 1 ( dC^\ T d 2 S a 



dlnT J v VPf V dlnV Jt Pv V dV ) T (Sy dVdT 
The Maxwell relation (dS cx /dV) T = (d(W/V)/dT) v = /3^ x allows us to rewrite this as 



(B6) 



<91nC^ x \ T fd(3 v < \ /dln/3' 



dlnT J v (3™ \ dT J v \ dlnT , , 
Because dlncy = dlnC™ at constant volume this implies that 



(B7) 



dln/3^ x \ (d\n.p 



Sinewy v \dlnCy* 



OCX 



= 1 . (B8) 



Thus at constant volume (3y is proportional to c™, i.e., 7 is constant on isochores. 

APPENDIX C: RELATING THE CORRELATION COEFFICIENT OF FIG. 1 TO THE WU 
CORRELATION COEFFICIENT R OF EQ. Q 

In order to distinguish from the WU correlation coefficient R of Eq. for the "direct isomorph check" of Fig. 
([IJ we denote by Rdj the correlation coefficient between potential energies of microscopic configurations with same 
reduced coordinates. Writing for a strongly correlating liquid for each microstate AW = ^AU + e, where e is an 
"error" term uncorrelated with AU, implies by squaring and averaging 

((AW) 2 )= 7 2 ((AC/) 2 ) + ( £ 2 ). (CI) 
Multiplying AW = 7 AU + e by AU, averaging and squaring gives 

(AWAU) 2 = 7 2 ((AJ7) 2 ) 2 . (C2) 
The left hand side of this is expressed in terms of the WU correlation coefficient R , giving 

((AW) 2 ) = ((AU) 2 ) / R 2 . (C3) 



Eliminating ((AW) 2 ) between this and Eq. (CI I gives a relation between the correlation coefficient and the variance 
of the error term, 

■i- - 1 = (£2) fC4) 
R 2 1 2 ((AU) 2 ) ' 1 ' 

Next consider an infinitesimal rescaling of all microscopic configurations to a new density, calling the old and new 
potential energies of corresponding microscopic configurations U^ and U^ 2 \ respectively. If the relative density 



change is small, Eq. ( Al I minus its average implies for each microscopic configuration 

AU {2) = AU {1) + (dlnp)AW il \ (C5) 

which, using AW = jAU + e, becomes 

AU {2) = (1 + 1 dlnp)AU {1) + (dlnp)e. (C6) 

We can now calculate the correlation coefficient Rdi between the old and new potential energies of microscopic 
configurations with the same reduced coordinates. In fact, since this equation is identical in structure to AW = 



7AC/ + e, we can use the result of Eq. (C4), replacing 7 with 1 + 7<ilnp and e with (dlup)e: 
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R 2 



(e 2 )(dlnp) 2 
(l + 7 dlnp) 2 ((A[/) 2 ) 



(C7) 



Using Eq. (C4| to eliminate (e } / ((AU) ) finally implies to lowest order in dlnp: 



r2 di 



— 1 = (7 din pY 



1 



- 1 



(C8) 



Equation (C8 1 applies to any liquid. If the density change goes to zero (dlnp — > 0) and/or the liquid becomes 100% 
correlating (R — > 1), there is perfect correlation between the new and old potential energies (Rdi — ► 1), as expected. 
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